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ABSTRACT 


The Thermal and Electrical Conductivity Probe (TECP) is a component of the 
Microscopy, Electrochemistry, and Conductivity Analyzer (MECA) payload on the 
Phoenix Lander. TECP will measure the temperature, thermal conductivity and 
volumetric heat capacity of the regolith. It will also detect and quantity the population of 
mobile H 2 O molecules in the regolith, if any, throughout the polar summer, by measuring 
the electrical conductivity of the regolith, as well as the dielectric permittivity. In the 
vapor phase, TECP is capable of measuring the atmospheric H 2 O vapor abundance, as 
well as augment the wind velocity measurements from the meteorology instrumentation. 
TECP is mounted near the end of the 2.3 m Robotic Arm, and can be placed either in the 
regolith material or held aloft in the atmosphere. 

This paper describes the development and calibration of the TECP. In addition, 
substantial characterization of the instrument has been conducted to identify behavioral 
characteristics that might affect landed surface operations. The greatest potential issue 
identified in characterization tests is the extraordinary sensitivity of the TECP to 
placement. Small gaps alter the contact between the TECP and regolith, complicating 
data interpretation. Testing with the Phoenix Robotic Arm identified mitigation 
techniques that will be implemented during flight. A flight model of the instrument was 
also field tested in the Antarctic Dry Valleys during the 2007-2008 International Polar 
year. 
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I Introduction 


Phoenix will be the first Martian lander to operate at polar latitudes, affording a unique 
opportunity to study the current climate, the role of surface-atmosphere exchange, and the 
polar climate record in the structure and distribution of ground ice. It will make the 
5 first-ever measurements to constrain volatile and energy fluxes in a region that is key to 
Martian climatic evolution, and the first physical and chemical analyses of Martian 
ground ice. By combining observations of present-day processes, and the geologic record 
of past processes, Phoenix will provide fundamental insights into Martian history, and its 
astrobiological potential. 

10 The Thermal and Electrical Conductivity Probe (TECP), a component of the Microscopy, 
Electrochemistry, and Conductivity Analyzer (MECA) on the Phoenix Lander, will 
address these topics. Based on the small, dual-probe sensors that are routinely used to 
monitor soil thermal properties and water content, TECP performs 6 distinct 
measurements: dielectric permittivity, electrical conductivity, temperature, thermal 
1 5 conductivity, volumetric heat capacity, and relative humidity. In addition, application of 

the thermal properties protocol when the TECP is held aloft by the Robotic Arm (RA) 
allows constraints to be placed on the average wind speed. 

This paper describes the TECP, its measurements, its design and implementation, its 
calibration, characterization and field testing, as well as some of the anticipated science 
20 results. The importance of the TECP dataset is enhanced by the uniqueness of the 
Phoenix landing site in the Martian polar regions. 



A. The High Martian Latitudes 

The Phoenix lander touched down on May 25, 2008 at 68.16°N 233.35°E, (Figure 1 and 
Arvidson et al., This Issue), an area that is fairly typical of the Vastitas Borealis plains. 
25 Investigation of a near-polar landing site gives us the opportunity to characterize a part of 
Mars that is more dynamic, and arguably, over the past several 10 7 years at least, more 
globally significant, than previous landing sites. Many components of the Martian 
climate are effectively controlled, or play out most dramatically , at these latitudes. 

H 2 0 Reservoir 

30 The high latitude surface is underlain by probably the largest H 2 O reservoir on the planet. 
Substantial permafrost deposits were first predicted on energy balance considerations 
(Leighton and Murray, 1966). Post-Viking, there was supporting geomorphic evidence 
for wide spread ground ice at high latitudes (Rossbacher and Judson, 1981; Squyres 
1986). The presence of high-latitude ice deposits were confirmed by Mars Odyssey 
35 (Boynton et al., 2002). 

Among the criteria for choosing the Phoenix site was that the modeled Odyssey GRS data 
suggests that the dry overburden thins somewhat at these longitudes; access to both the 
ground ice and its dry overburden were considered fundamental objectives. 

Seasonal Cycle 

40 The Phoenix landing site is also part of the circumpolar terrain covered annually by the 
seasonal CO 2 cap. Within the seasonal cap perimeter, the surface energy balance controls 
the condensation of H 2 O and CO 2 . Cap condensation and sublimation are the dominant 
climatic processes in the current epoch. Phoenix is well placed to observe and begin 
characterization of cap condensation processes. 
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45 Once the seasonal cap has retreated poleward of the Phoenix site, the high-latitude 
regolith is likely a source of H 2 O to the atmosphere, which along with cap, buffers 
planetary atmospheric H 2 O (Jakosky, 1983). Depending on the thermal and diffusive 
properties of the regolith, H 2 O vapor exchange with the ground ice could have a major 
impact on the annual atmospheric water budget (Schorghofer and Aharonson, 2005). 
50 Results from the Planetary Fourier Specrrometer (PFS) suggest the possibility that 
atmospheric H 2 O is preferentially concentrated in the lowest few kilometers of the 
surface at high northern latitudes during summer (Tschimmel et al., 2008). If true, this 
would affect the upper boundary condition in atmospheric H 2 O, and hence subsurface ice 
stability models. 

55 Obliquity Effects 

In addition to the important role the polar regolith plays in determining seasonal cycles, 
the high latitudes experience the most extreme variations over orbital oscillations. In 
particular, variations in orbital obliquity have a pronounced effect on polar temperatures. 
At low obliquity, variations in polar insolation over the year are minimized. It was 
60 formerly believed that semi-permanent CO 2 caps would form at both poles at minimum 
obliquity (Toon et al. 1980, Pollack and Toon, 1982, Fanale et al, 1982), cold-trapping 
regolith ice away. More recent modeling however, calls into question the nature of the 
cap at low obliquity (Haberle et al., 2003; Armstrong et al., 2004). 

When obliquities are high, the polar regions experience extreme seasonal insolation 
65 variations. The north polar H 2 O cap becomes unstable at some obliquity between 35° and 
45° (Mischna et al., 2003). The stable location for ice at obliquities > 45° is the low- 
latitude highlands, as first suggested by Jakosky and Carr, (1985), but more recently 
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supported by Head et al. (2003). 

Because it is at the poles that obliquity variations have the most pronounced effect, the 
70 subsurface at the Phoenix site may provide a record of the processes and mechanisms that 
drive long term climate change. 

Periglacial Processes 

There are other similarities between the Martian high latitude surface, and some of 
Earth’s periglacial terrain. Understanding their true origin and nature is essential if we 
75 are to read the geologic record, both at the surface, and in the subsurface. 

Polygonal Ground 

Much of the Northern Plains, including the Phoenix landing site, are covered with 
patterned ground that is remarkably similar in appearance to examples that are commonly 
seen at terrestrial polar latitudes. (Figure 1; Arvidson et al., This Issue ; McEwen et al. 
80 2007) 

Polygon formation begins with formation of a network of vertically-propagating cracks 
due to the thermal contraction of ground ice during autumn and winter. Mellon (1997) 
modeled the process for the Martian climate, and found the yield strength of frozen 
ground is exceeded during autumnal cooling, leading to mechanical failure and crack 
85 formation (Lachenbruch, 1962; French, 1996). The cracks then accumulate either clastic 
particles, or, if it is available, liquid water. It is also possible to condense additional H 2 O 
vapor in these cracks, a process that could likewise contribute to the accumulation of 
excess ground ice (Fisher, 2005). 

At current pressures and temperatures, liquid water is not stable on Mars long enough to 
90 fill the cracks, but past conditions may have permitted this process (Costard et al., 2002). 
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Whether these features are relict of freeze-thaw cycles, due solely to thermal contraction, 
or are the product of other processes is a key question that Phoenix must resolve. 

Excess Ground Ice 

The Gamma Ray Spectrometer on Mars Odyssey shows the presence of an ice-rich 
95 layers with unexplained abundances of more than 60% of ice by volume (Boynton et al., 
2002; Feldman et al., 2002). Buried H 2 O deposits, such as snow, layered with dust 
(Mischna et al., 2003) are a possible explanation for the large amount of interstitial ice, 
but periglacial ffeeze-thaw cycles could produce the same result.. 

During freezing events, as the freezing front penetrates downward into the surface, 
100 unfrozen, interfacial water is removed to the ice phase, creating a hydraulic gradient 
along which additional water migrates. As long as the total (sensible + latent) heat flux 
carried by unfrozen water from below can balance the diffusion of heat upwards to the 
surface, the freezing front will remain stationary, and ice lenses will develop in freezing 
soil. Once the regolith beneath the freezing front is depleted of H 2 O, the sensible + latent 
105 heat flux to the freezing front is no longer sufficient to balance conduction, and the 
freezing front again begins to propagate to depth. This continues until another wet zone 
is reached, and another ice lens begins to form. Terrestrial periglacial terrain often 
exhibits numerous ice lenses in the shallow subsurface. 

Habitability 

110 The thermal physics of Martian permafrost are also of interest because Earth’s high- 
latitude permafrost, which bears many morphological similarities to the Phoenix landing 
site, is a rich habitat for psychrophiles. 

Rivkina et al., (2000) argue, based on growth metabolism in the laboratory, that if the 
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115 


temperature is stable, as in deep permafrost, the minimum temperature for growth is 
perhaps 253 K. Growth metabolism is easily detectable at 263 K (Finegold, 1996). 

Price and Sowers (2004) reviewed the literature on psychrophilic life, and could find no 
detectable lower temperature limit for metabolism. They concluded that the extremely 
low energy costs of survival energy enable microbial communities in extreme 
environments to survive indefinitely. They also reported evidence for survival-level 
120 metabolism all the way down to 23 3K, meaning that much of the Martian regolith and 
ground ice qualifies as “survivable” if only thermal limits are considered. 

Jakosky et al (2003), examined the minimum temperatures at which interfacial unfrozen 
water can exist, the minimum subfreezing temperatures at which terrestrial organisms can 
grow and multiply, and the maximum temperatures that can occur in Martian high- 
125 latitude and polar regions, to see if there is overlap. They concluded that high-latitude 
Martian ground ice was just at the limit of being habitable. 

The most probable lifetime limitation is exposure to galactic cosmic rays (Kminek et al., 
2003); an accumulation of damage eventually overwhelms the cell’s repair mechanisms. 
Kminek et al (2003) predicted that the time required to reduce a viable spore population 
130 to 10' 6 would be < 6xl0 5 years in the top meter of the regolith, and < 7xl0 4 years in the 
otherwise-favorable upper 10 cm. 

Zent (2008) used a time-resolved thermal model of the subsurface, tuned to the Phoenix 
landing site, along with the 10-My orbital calculations of Laskar et al. (2004), to calculate 
several possible 10 My thermal histories of Phoenix’s target ice. In all cases, there were 
135 many clement periods when growth, even rapid growth, was possible. Also however, in 
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every case, the maximum dormancy period that must be survived was at least an order of 
magnitude greater than the ~7xl0 4 years quoted by Kminek et al. (2003) 

B. Science and Measurement Objectives 

The science objectives of the Phoenix Mission are: 

140 1 . To study the history of Martian water in all its phases. 

2. To search for evidence of habitable zone and assess the biological potential of the ice- 
soil boundary. 

The polar regolith holds a substantial fraction of the remaining Martian H 2 O inventory, 
and is a dynamic reservoir for H 2 O over both seasonal and astronomically-driven climate 
145 cycles. The abundance, distribution, and structure of polar ground ice may also have 
preserved a geologic record of the recent climate and subsurface conditions. Finally, as 
we have already discussed, polar ground ice may be just at the margin of habitability. 

The TECP was included on Phoenix specifically to determine the regolith thermophysical 
properties that control energy fluxes between the atmosphere and subsurface, to 
1 50 characterize the processes that control the distribution and exchange of H 2 O between the 
atmosphere and subsurface in the current climate, and to document the occurrence of any 
unfrozen water, both as thin films in the regolith, and as atmospheric vapor. 

1 . Science Objectives 

a. Energy Fluxes 

155 Regolith thermal properties, particularly the thermal conductivity (k) and the volumetric 

heat capacity (Cy) are key, and currently unknown, parameters in any model that purports 
to address the distribution of Martian ground ice, the climatic cycles associated with 
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variations in orbital parameters, or the annual cycles of H 2 O and CO 2 . These properties 
determine how effectively heat conducts from the illuminated surface to depth, and fix 
1 60 the thermal profile from the surface down at least as far as the ice table. 

Based upon these measurements, estimates of the thermal inertia of the surface can be 
made, and compared with thermal inertias that are derived from measurement of the 
surface temperature throughout a thermal cycle, whether by TECP itself, or with thermal 
inertias derived from orbital data (e.g. Mellon et al., 2000). 

165 b. Regolith H 2 0 

Another set of objectives involves understanding the abundance, state, and mobility of 
H 2 O in the regolith. H 2 O redistributes itself continuously through the regolith, 
continually accumulating in the lowest-energy reservoir that is thermodynamically 
available. Because soil and ice co-exist in the regolith, it is necessarily true that there is 
170 a finite and variable population of un-frozen H 2 O molecules; if sufficiently thick films of 
unfrozen water occur, surface diffusion in the film can outstrip diffusion in the vapor 
phase. These mobile H 2 O molecules are most easily detected in measurements of soil 
electrical properties, which are otherwise dominated by highly insulating silicates and 
ice. 

175 Because TECP can measure both unfrozen water and temperature over a diurnal cycle, it 
may be possible to characterize both the soil water characteristic curve (the relationship 
between soil water content and soil water pressure potential) as well as the hysteresis 
curve of H 2 O adsorption in the regolith materials. These are key parameters in 
determining the flow and transport of unfrozen water in the regolith. 
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180 c. Atmosphere 


The atmosphere serves as the conduit through which H 2 O is distributed around the planet, 
and as the upper boundary condition controlling the state and distribution of H 2 O 
throughout the regolith. Exchange between the atmosphere and regolith is controlled 
substantially by the absolute humidity of the atmosphere, and somewhat to a lesser extent 
1 85 by the wind fields. Therefore, TECP will also measure atmospheric relative humidity, 
and, by applying the same protocol used to determine regolith thermal properties, can 
constrain average wind speeds averaged over the duration of the measurement. 


2. Measurement Objectives 


190 


195 


200 


The measurement objectives for TECP are: 

• Measure the temperature of the Martian regolith, both at the surface, and within 
the excavation trench, as a function of Local Solar Time (LST), and the sub-solar 
angle (L s ) which is used to express Martian seasons. If different types of 
materials can be accessed, these measurements will be repeated over the same 
timescales for each material. Based upon surface temperature observations 
(Kieffer and Titus, 2001), measurement requirements on the temperature range 
from 165K to 323 K, with 2K precision, and an accuracy of ±10%. 

• Measure the thermal conductivity and heat capacity of accessible Martian regolith 
materials. The TECP thermal conductivity requirement is to measure soil thermal 
conductivity (k) with 10% accuracy over the range 0.03 - 2.5 [W m' 1 K' 1 ]. A 

typical value for thermal conductivity in frozen soils or solid ice is k= 2.0 ± 1.0 
[W m' 1 K' 1 ] at T= 223 K. This value is from many earth observations over many 


11 



205 


210 


215 


220 


types of soil, including pure ice. We concluded that we would not require TECP 
to measure pathologically low thermal conductivity, since the landing event and 
the measurement attempt itself would probably destroy any such fluffy surface; 
simply confirming an off-scale low K would be useful. Note that for such a 

surface the thermal inertia, 7= <Jkcp = (0.03*800 [J/kg K]*800 [kg/m 3 ] = 138, 
which is extremely low. The upper limit on the thermal conductivity requirement 
was a result of heritage design constraints. Thermal conductivities above 2.5 [W 
m' 1 K' 1 ] result in very small temperature increases at the heated needle, which are 
difficult to resolve with the 12 bit ADC. 

• The volumetric heat capacity (Cy) of a particulate mixture is the heat capacity of 
each component, weighted by its volumetric relative abundance; or approximately 
the specific heat of the regolith times its density. The requirement for the 
instrument is to measure the volumetric heat capacity of the regolith with 10% 
accuracy over range 0.4 to 4 [MJ/m K], which approximately spans the range 
from dry frozen ground to nearly the heat capacity of liquid H 2 O (4.19 [MJ/m 
K]) (Hinkel et aL, 2001; Kane et al., 2001). 

• Measure the electrical conductivity (inverse of the electrical resistivity) of the 
regolith materials, particularly as a function of measured regolith temperature. 
The requirements for electrical conductivity are to measure soil and ice electrical 
conductivity from 10 7 [nS/cm] to 1 [nS/cm] with 10% accuracy. Conductivity, (a) 

is a property of a material (e.g., copper), as opposed to conductance, which is a 
property of a specific conductor (e.g. a wire made of copper). Conductivity has 
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225 


230 


235 


240 


245 


units of Siemens per meter [S/m]; conductivity will typically be referred to in 
terms of [nS/cm] 


a 



w 

m 2 J 


= 10 7 


nS 

cm 


[ 1 ] 


De-ionized water has a conductivity of 5.5xl0" 6 [S/m] (55 [nS/cm]); Siberian 
permafrost has been measured with electrical resistivity, R, as high as 1.4xl0 4 [Q- 

m] (a = 715 [nS/cm]) at depth > 2.2 m just before warming begins (Olhoeft, 
1978; Pozdnyakov and Chan, 1979). 

• Measure the bulk dielectric permittivity of the regolith materials over the range of 
1 to 20, with resolution of 0.005, as a function of regolith temperature. The 
dielectric permittivity of any gas is ~1, and £ for water is 78.54. We chose not to 

measure the full range because high sensitivity (hence high resolution) is desired 
at the low end of the range to maximize sensitivity to small changes in the 
abundance of unfrozen water. The upper limit of 20 is imposed because MECA 
data words are only 12 bits, so achieving sensitivity of 0.005 at the low end of the 
range limited the upper limit of the dynamic range to 20. 

The relationship between the water content of a soil and dielectric permittivity is 
complex and continues to be studied (Topp et al. 1980; Jones and Or, 2003; 
Lebron et al, 2004; Sweeney et al. 2007). The overall dielectric is the sum of the 
dielectrics of air, mineral, ice and water. In a non-frozen system, we typically 
assume unfrozen water is the only constituent changing in the medium; all other 
dielectrics stay the same. The change in dielectric permittivity can then be 
interpreted solely in terms of the change in unfrozen water. However, T changes 
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will affect permittivity as well; therefore all £ measurements are carried out with 


simultaneous T measurements of the medium. 

• Measure the H 2 O vapor density of the Martian atmosphere as a function of LST, 
L s , and height above the surface, including measurements made within the trench 

250 and particularly just above the ice layer if an ice layer is detected. TECP also 

includes a GE Panametrics MiniCap 2 polymer relative humidity (Rh) sensor, a 
commercial capacitive polymer-based sensor whose performance at temperatures 
and frost-points below 233 K was not well known initially (see the Calibration 
Discussion, below). The requirement is to measure the H 2 O vapor pressure from 0 
255 - 0.5 kPa over the temperature range of 1 95 K - 270 K, with accuracy of ±10% of 

reading. Because the sensor is sensitive to Rh, the absolute partial pressure can 
only be determined with a measurement of air temperature. Strategies for 
determining absolute H 2 O partial pressure are discussed below. An additional 
consequence of this approach is that the resolution in absolute pressure is a 
260 function of the sensor temperature, which is related to, but not identical to, the air 

temperature. 

• Measure, or characterize the wind speed as a function of LST, L s , and height, 
including within the trench. 

II Technical Approach 

265 In order to meet the measurement requirements, TECP was developed as a single 
electronics box, (Figure 2) integrated with four needle-like “probes”, each of which has a 
role as a component in a unique set of circuits (Figure 3). 
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The TECP needles are fitted into a plastic housing, which in turn is attached to the 
instrument electronics box. The plastic interface is machined from polyethylethylketone 
270 (PEEK), which has low thermal and electrical conductivity (k = 0.249 W m' 1 K' 1 ; Cy = 

2. Ml kJ kg' 1 K' 1 ), and provides insulation between the thermally and electrically 
activated needles and the electronics box. The aluminum electronics housing is anodized 
to increase emissivity, enhancing thermal radiation, and minimizing heating when TECP 
is operating in direct sunlight. A port in the side of the instrument housing allows the 
275 relative humidity sensor to access the Martian atmosphere. The sensor is shielded behind 
a Teflon membrane that is porous to atmospheric H 2 O but prevents dust from fouling the 
humidity sensor. 

TECP has a single rotational degree of freedom relative to the RA, (Figure 4), which in 
turn has 3 degrees of freedom relative to the lander and the environment. 

280 A. Temperature 

The temperature of three of the needles (1, 2, and 4) is measured as the difference 
between the sense and reference junctions of a Type E (chromel-constantan) 
thermocouple. The thermocouple was fabricated by soldering the ends of a 7.5 cm piece 
of constantan thermocouple wire to the ends of two 75 cm pieces of chromel 
285 thermocouple with Kester Sn63Pb37 rosin core solder and Indium Flux #1. The reference 
junction of the thermocouple is on the TECP analog electronics board, adjacent to a 
National Semiconductor LM234DT temperature-sensing current source. The temperature 
derived from this element is referred to as the board temperature (I*); the output current 
is directly proportional to the absolute temperature, and extrapolates to zero at 0 K. 
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290 B. Thermal Properties 

The in situ thermal properties of the Martian regolith are determined via the transient 
heated needle technique (de Vries, 1952). In this approach, a heat pulse is applied via a 
heater that approximates a line source, and the thermal response of the medium is 
measured at some radial distance from the heat source. 

295 In addition to thermocouples, needles 1, 2 and 4 are also equipped with resistance 
heaters. The heater is Evanohm, an alloy used extensively for precision wound resistors. 
Its composition is nominally 0.75 Ni, 0.20 Cr, 0.025 Al, and 0.025 Cu. The Evanohm is 
dipped in lacquer because the thermal grease used to ensure good thermal contact to the 
needles has a metallic component, which would cause the heater to short to the needles. 
300 In typical operations, needle 1 is heated, and the temperature change of both needle 1 and 
needle 2 is monitored., The temperature change of needle i, between time 0 and time t is 
determined from 

AT t>t =(T i>t -T 4>t )-(T ifi + T 4fi ) [2] 

Needle four (farthest from the heated needle) is used as a temperature reference, as it is 
305 unaffected by the heat pulse from the needle for heat pulses of normal duration. This 
method also has the advantage of removing the effects of thermal drift of the bulk 
material or the TECP electronics board during the measurement. 

Carslaw and Jaeger, (1959) modeled the temperature surrounding an infinite line heat 
source with constant heat output and zero mass, in an infinite medium. When a constant 
310 amount of heat is emitted from unit length of a zero-mass heater over a period of time, 
the temperature response at radial distance r, is 


16 



0<t<ti [3] 


A T tl =--^—Ei 
■ 4tuk 


^-r 2 C } ' 

4ict 


where q is the rate of heat dissipation [W m" 1 ], t\ is the heating time, and Ei is the 
exponential integral. The temperature change after the heat is turned off is given by: 


315 


A T tyi = - 


4 K K 


-Ei 


f -r 2 C, ' 

v 4Kl y 


+ Ei 


r -r 2 Cy 
4K{t-t\) 


t>h [4] 


Material thermal properties are determined by fitting the time series temperature data 
during heating to Eq. (3), and during cooling to Eq. (4); note that the absolute 
temperature T doesn’t enter into the calculation, only AT.. Thermal conductivity and heat 

capacity can be obtained from the temperature of the heated needle (single needle), with r 
320 taken as the radius of the needle and independently by fitting equations 3 and 4 to the 
temperatures measured at the adjacent needle, 7 mm from the heated needle. 

For the fitting procedure, a multi-stage Monte Carlo method described by Conley (1985) 
is used to obtain values close to the global minimum, and then a Marquardt method is 
used to find the final values. The exponential integral function was approximated using 
325 formulae from Abramowitz and Stegun (1972). 

The outputs of the solver function are k ^ and C/i, which are the thermal conductivity 

and volumetric heat capacity measured at the heated needle (needle 1), and k / 2 and 

C/2, which are thermal conductivity and volumetric heat capacity measured at the 

adjacent needle (needle 2). The values for k and C/ found via fitting to Equations 3 

330 and 4 are inherently inaccurate due to inevitable discrepancies between the TECP and the 
mathematically idealized model. However, the true thermal properties can be determined 
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if functional calibration is performed against materials of known thermal properties as 
described below. 


C. Dielectric 

335 To measure the dielectric permittivity of the regolith, two of the needles (3 and 4) are 
used as the plates of a capacitor; with the regolith serving as the dielectric. The 
measurement involves putting a microwave frequency (~6.25 MHz) voltage square wave 
on the capacitor, and measuring the time history of the charge on the plates. 

The time, t, required to charge a capacitor from an initial voltage (Vi) to a final voltage 
340 (V/), with an applied voltage (V) is 


t = -/?Cln 


f V f - V 


V t -V 


[5] 


where R is the series resistance and C is the capacitance. The capacitance, in turn, is a 
function of the dielectric permittivity of the medium between the plates (or in this case, 
the regolith between the needles) 

345 C = eZc [6] 


where Lc is a geometric factor that encompasses the area of the capacitive plates and their 
separation. Lc has units of length [m]. The relationship between the time history of the 
plate charge and the dielectric permittivity is then 


1 

e 


( 


-RL cln 


V f -V 

v-v 


y 

). 


[7] 


350 Because we are measuring the effective permittivity of a mixture of phases, rather than a 
pure phase, we refer to the measured quantity as the bulk dielectric permittivity (£b). 
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D. Electrical Conductivity 

The electrical conductivity (a) of the soil is measured between needles 1 and 2. The 
regolith conductance ( Gr ), the inverse of the regolith resistance (Rr), is measured via a 
355 simple voltage divider (cf. Figure 3). Three fixed resistors (Rf = lkQ; 100 kD; 3MQ) are 

used to allow the 12-bit data word to cover the entire dynamic range specified in the 
requirements with sufficient precision. This results in three distinct measurement 
channels, each with its own effective range, which are measured sequentially (Table 4). 
The input voltage (Vi) is 2.5V, and the resistance of the regolith is found from the output 
360 voltage ( V 0 ) via: 


f \ 


Rr = 


V 0 


Rf 


V 



[98] 


The conductivity of the regolith is found by multiplying conductance (1 /Rr) by the probe 
constant ( Cp ), which in effect normalizes for the geometry of the electric field that 
permeates the regolith between the needles. The determination of the probe constant is 
365 described below. 


E. H 2 0 Vapor 

The TECP H 2 O measurement is based on a GE Panametrics MiniCap 2 polymer relative 
humidity sensor. It is a capacitance-based instrument that measures the permittivity of a 
polymer film which adsorbs H 2 O. The circuit in the TECP is read out in a manner similar 
370 to the capacitance measurement already described for the soil measurement. 

Solid-state humidity sensors of this type are sensitive to relative humidity, almost 
independent of temperature. Specifically, the Rh is measured relative to the equilibrium 
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pressure of liquid H 2 O, rather than ice, regardless of temperature. At T < 273 K, the 
saturation vapor pressure of supercooled water is used. 

375 The physical basis for this behavior is likely that atmospheric H 2 O forms an adsorbed 
film on the sensing polymer in equilibrium with the atmosphere, and it is the thickness of 
this adsorbate that is actually sensed as a change in permittivity It can be shown that, for 
many adsorbents, the mass (or thickness) of this adsorbed layer is nearly the same for a 
given R h , regardless of temperature (Anderson, 1995). Above some limiting film 
380 thickness, corresponding to Rh = 100%, the thickness of the film increases without limit, 
and bulk H 2 O condenses. 

Unfortunately, there is a low-temperature limit in humidity measurements of this type. 
To substantially affect the bulk permittivity of the sensing film, the H 2 O molecules must 
be free to align with the applied electric field; that is they must remain “unfrozen”. 
385 Should some of the water adsorbed in the sensing film undergo spontaneous 
crystallization to ice, the device would instantly lose sensitivity to that portion of H 2 O, 
resulting in a discontinuity in sensor output. 

One may think of the interfacial, unfrozen water in the sensor as being in dynamic 
equilibrium with the competing forces of the H 2 O hydrogen bonds in ice, and the van der 
390 Walls bonds of adsorbing polymer. As T decreases below 273 K, the hydrogen bonds 
between H 2 O molecules increasingly overcome the sum of the H 2 O kinetic energy + the 
energy of the lUO-polymer bond, and the maximum possible thickness of the unfrozen 
water layer decreases. At whatever temperature this ever-decreasing maximum reaches 
the film thickness corresponding to Rh — 100%, ice will condense, and the instrument 
395 will begin to lose sensitivity, beginning at the highest humidity. Experience suggests that 
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in the TECP, this lower limit is reached ~ 205 K. Additional characterization work, and 
examination of the flight data will determine this limit more precisely. 

In TECP, the humidity sensor is located adjacent to the board temperature sensor on the 
analog electronics board, and Tj is used a a measure of the Rh sensor. Due to power 
400 dissipation on the electronics board, is typically slightly higher than the ambient air 
temperature during operation. This works to our advantage, in that it maintains the 
sensor above the frost point when the atmospheric Rh is high. However, it also means the 
Rh of the sensor is lower than that of the atmosphere. In practice however, we can use 
the measured Rh along with Tb to calculate the absolute humidity of the atmosphere. If 
405 an estimate of atmospheric Rh is desired, it can be calculated from absolute humidity 
along with an independent atmospheric temperature measurement, such as those made by 
the meteorology payload. 

Ill Calibration 

A. Temperature 

410 1 . Board Temperature Sensor 

The board temperature sensors were calibrated to NIST-traceable temperature standards 
in a two part process. Opening the flight units and attaching temperature sensors directly 
to the boards was not an acceptable option, so an Engineering Model (EM) unit was used 
to characterize the temperature difference between 7), and the temperature at the outside 
415 of the TECP enclosure (on the aluminum cover directly above the board temperature 
sensor). NIST-traceable temperature references were then attached to the flight TECP 
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units above the board temperature sensors and measured while the TECP was brought to 
temperatures over its full expected temperature range. The correction determined on the 
EM unit was then applied and the board temperature calibration function was determined. 
420 The temperature difference between the PCB and the analog board cover varied smoothly 
as a function of temperature according to 

A7 7 = — 4.39 17x1 0 -6 /” 2 + 2.8044x10^7 + 0.64736 [9] 

with R = 0.996. This polynomial was applied to the calibration data. 

Diode type, NIST-traceable temperature sensors (Lake Shore Cryotronics model DI-470) 
425 were attached with kapton tape to the TECP aluminum cover directly over the board 
temperature sensors on the TECP flight units. The TECP units were placed in an 
aluminum thermal enclosure in a Delta Design environmental test chamber. The 
chamber was used to bring the thermal enclosure and TECP to equilibrium at eight 
temperatures: 160 K, 180 K, 200 K, 220 K, 240 K, 260 K, 280 K, and 323 K. During 
430 post-processing, the difference between T\, and the analog cover temperature (Eq. 9) was 
applied to the temperatures measured by the diode sensors. Figure 5 shows the calibration 
function for on the TECP flight unit. 

2. Thermocouple Calibration 

Thermocouple calibration involved determining the Seebeck coefficient (Cs) vs. 
435 temperature for the thermocouple alloys used in the TECP needles. This allows accurate 
determination of the temperature difference between the TECP needles and 7*. The 
Seebeck coefficient measured here was from the same spool used to fabricate the flight 
model thermocouples. 
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Two cylindrical pieces of copper 1.3 cm long x 4.6 cm diameter were bolted 6.5 cm apart 
440 on a piece of acrylic, which was set on top of a block of open cell foam in a Delta Design 
temperature chamber. The two blocks each had a 3.2 mm diameter hole bored to their 
center from the rounded surface which housed pre-calibrated temperature sensors and 
thermocouple junctions. The thermocouple junctions and NIST traceable temperature 
sensors were co-located and potted in the borehole in the copper blocks with Arctic Silver 
445 high thermal conductivity grease. 

In one of the blocks, a 0.4 cm hole was bored completely through the block to facilitate a 
180 Q electrical heating element, allowing the block temperature to be controlled 

independently, and thus achieve a temperature difference between the two thermocouple 
junctions. (Figure 6). 

450 The Delta Design temperature chamber was used to control the temperature of the 
atmosphere around the thermocouples. Seebeck coefficients were determined for the 
thermocouple at approximately 30 K increments from 150 to 330 K. The blocks were 
allowed to equilibrate at each ambient temperature setting for at least 30 minutes. At 
each ambient temperature, a voltage between 0 and 32 V was applied to the heater in the 
455 copper block producing a temperature difference (AT) between the two blocks. At least 
60 seconds were allowed to pass between the final voltage adjustment and the recording 
of data to ensure that the system was in thermal equilibrium. 

Figure 7 shows the Seebeck coefficient vs. temperature data collected here, which closely 
follow a third order polynomial function. The Seebeck coefficient function shown here is 
460 used to calculate the temperature difference between the needles and Tj, on the flight unit. 
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B. Thermal Properties 

The TECP thermal properties measurement functions were calibrated at the system level. 
Measurements in materials of known thermal properties were used to construct a 
calibration function for the thermal conductivity measurement and volumetric heat 
465 capacity measurement for TECP in a manner described in Cobos et al., (2006) and 
summarized here. 

TECP thermal properties time series data sets were collected in five materials with 
known thermal properties (Table 1). The Macor ceramic, high density polyethylene 
(HDPE), and glycerol were obtained from commercial vendors. The expanded 
470 polystyrene (EPS) was obtained from a larger block of certified material used as a 
thermal conductivity calibration standard at Decagon Devices. The water was stabilized 
with agar powder at a mixing ratio of 5 [g/1] to prevent convective heat loss from the 
heated needle. 

In the case of the stabilized water and glycerol, thermal properties values were taken 
475 from published reference literature. The thermal properties of the Macor and HDPE 
samples were determined with a commercial thermal properties meter (model KD2 Pro, 
Decagon Devices), which was calibrated immediately prior to testing. The EPS material 
is a commercially available, NIST traceable thermal conductivity standard material 
obtained from LaserComp, Inc. The volumetric specific heat of the foam standard was 
480 calculated calorimetrically from published values of the specific heat and density of 
polystyrene and the measured density of the EPS. 

All measurements were taken in a drybox using a vertical stage to insert the TECP 
needles into the standard materials. An N 2 purge was used during measurement of the 
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485 


glycerol, which is hygroscopic. The stabilized water and glycerol were contained in 600 
ml beakers, into which the TECP was lowered until the metal sensing needles were fully 
immersed in the fluid, but the PEEK shoulders were not. The TECP sensing needles 
were inserted into the EPS block using the vertical stage, and were again inserted until 
the needles were fully encased, and the PEEK shoulders were not. The HDPE and Macor 
samples were machined to accommodate the four TECP needles. A total of six needle 
490 holes were machined into the HDPE and Macor blocks, with the center two holes 
purposely cut slightly under-sized to ensure that the cones of needles 1 and 2 made solid 
thermal contact with the material. This resulted in a small gap between the PEEK 
shoulder and the calibration block. Before inserting the TECP into the Macor and HDPE 
blocks, all of the holes were completely filled with DI water to promote thermal contact 
495 between the needles and the material. The excess water that was forced from the holes 
during the insertion of the TECP was removed from the surface of the material block 
with a pipette. After insertion into each material, the system was allowed to come into 
thermal equilibrium before data were collected. Additionally, between repeated 
measurements on the same sample, the system was also allowed to come back to thermal 
500 equilibrium after introduction of the preceding heat pulse. 

For each thermal properties measurement, the TECP control software was configured to 
record the output of the board temperature sensor, the three needle thermocouples, and 
the heater current sense resistor at 1 Hz for 160 s. A 60 second heat period was applied 
to needle 1 after an initial 10 second buffer period with no heat. After cessation of the 
505 heat pulse, 90 s of cooling data were recorded. Each thermal properties measurement 
was repeated at least twice in each material. For the glycerol, the needles were removed 
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and re-inserted into the same sample between measurements. For the stabilized water, 
the needles were removed from the first sample and inserted into an identical, virgin 
sample between measurements. For the EPS, the needles were removed and inserted into 
510 a virgin section of the same block of the material. For the HDPE sample, the needles 
were removed from the sample, and the orientation of the sample was reversed so that 
needle 1 was in the hole formerly occupied by needle two. For the Macor sample, two 
measurements were made in the original orientation, and two in the reversed orientation 
for a total of four measurements. 

515 The measured thermal conductivity (k ot ) and heat capacity (C m ) are related to the values 
derived from the fitting process by the following calibration functions: 

K m = Q K/1 + C 2 Ky 2 + C 3 Cy! + C 4 Cf 2 + C 5 [10] 

and 

C m = Q K /1 + C 7 Ky 2 + CgCyi + C 9 Cf 2 + Qo [1 1] 

520 The calibration coefficients C 1 ....C 10 were determined during calibration by minimizing 
the error between actual values (k q and C a ), and the measured values (K m and C m ) 

determined in Eqs. 10 and 11. Values for Kyi, K \fi, Cj\ and Cfi, generated from the 

MSMCO-Marquardt fitting process for the 5 materials tested are shown in Table 2. The 
relationship between the measured values, derived from equations [10] and [11], and the 
525 actual values for the 5 known calibration materials results in the relations shown in 
Figures 8 and 9, which are the final conversion from the measured to actual k and Cy 
values. 
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C. Dielectric 


TECP dielectric measurements were taken in a series of four fluids with relative 
530 dielectric permittivities from 1 to 20 (Table 3). Each of the fluids was housed in a 30.4 
cm diameter by 17.5 cm depth glass bowl, filled to within 1 cm of the top. Earlier finite 
element analysis and experimental testing showed that >95% of the electromagnetic field 
generated by the TECP dielectric measurement was contained within this volume for all 
calibration fluids. All measurements were performed in a drybox under dry N 2 purge; the 
535 TECP was lowered into the fluids with a translation stage. The TECP units were inserted 
into the fluid until the needles were submerged. For the air, isopropyl, and 1-decanol 
samples, the fluids were allowed to equilibrate to the temperature of the drybox for at 
least 1 hour prior to data collection. The tert-butanol sample was held at a steady, above- 
ambient temperature with an electrical heating pad due to its high melting temperature 
540 (301 K). In each fluid, data were collected from the TECP unit for 25 minutes to assure 

thermal equilibrium between the needles and the fluid. Thermal drift was not identified 
in any of the samples after one minute of equilibration time, so the TECP output values 
presented are the averages of approximately the final 24 minutes of data collection. 

Prior to making TECP measurements in the fluids, six sub samples of each of the three 
545 measured alcohols were collected. The fluids were transferred to 100 ml, septa-equipped, 
airtight glass containers in the drybox under N 2 purge. The relative dielectric permittivity 
of the standards was independently determined at the National Institute of Standards and 
Technology’s Dielectric and Magnetic Material Laboratory. The NIST measurements 
were performed with two different parallel plate devices at the same measurement 
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550 frequency (6.25 MHz) and temperature as the TECP calibrations. Using the data in Table 
3, the TECP calibration function was determined (Figure 10). 

D. Electrical Conductivity 

The electrical conductivity calibration was a two-step process. First, the TECP 
conductivity function was characterized with a set of fixed resistors. Second, aqueous 
555 standard solutions were used to find the probe constant vs. resistance. By combining 
these two functions, the electrical conductivity of the bulk medium surrounding the 
needles can be accurately determined for each of the three resistance channels (Table 4). 
The TECP flight unit was placed in the Delta Design chamber. Copper blocks machined 
at the same cone angle as the TECP needles were used to make electrical contact between 
560 the needles and two clip-type leads that were routed through mouseholes to the exterior 
of the chamber. The resistance of the leads was measured at 0.06 Q, and has been 
neglected in the following analysis. The clip leads were used to sequentially connect the 
TECP needles across a series of 24 axial lead resistors ranging from 50 Q to 250 MQ. 
The resistance of each of the fixed resistors was measured just prior to each TECP 
565 measurement with a calibrated Agilent multimeter. To allow circuit settling, data were 
collected across each resistor for 120 seconds for each of the three measurement ranges. 
The average of the last 10 seconds of TECP output for each measurement range has been 
used in all subsequent analyses. This process was repeated at five different temperatures, 
spanning the range 160 K to 323 K. 

570 Measurements were then taken in a series of six aqueous standard solutions ranging from 
1 to 10,000 qS/cm. The standards were poured to 500 ml in a 600 ml Pyrex beaker. All 
experiments were performed in the dry box with the TECP lowered into the solutions 
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with the translation stage until the metal needles were fully submerged. The PEEK 
shoulders were not submerged in the solution. Measurements began with the 1 |_iS/cm 

575 sample (DDI water), and progressed to the highest salt concentration sample. The TECP 
needles were wiped dry and wiped with an isopropyl alcohol clean room wipe between 
samples. Before measuring with the TECP, the conductivity of each solution was 
measured with a calibrated Orion multi meter. In the special case of the 1 |_iS/cm 

solution, where a very small amount of salt contamination would affect the conductivity 
580 of the standard significantly, the TECP was first dipped into the sample and removed, 
then the conductivity was measured with the Orion meter, and then the TECP was placed 
back into the sample and a measurement was taken. The calibration again consisted of a 
120 second measurement in each of the three measurement channels. Again, the average 
of the last 10 seconds of TECP output for each measurement range has been used in all 
585 subsequent analyses. 

The data collected across the 24 resistors were used to construct a matrix of ADC output 
as a function of resistance and temperature for each channel. After each data set was 
abbreviated to its useful range (cf. Table 4) it was modeled as a 4 th order polynomial 
function of TECP output vs. In resistance (Table 5). 

590 Preliminary characterizations indicated that the TECP probe constant (Cp) varies as a 
function of the measured resistance at least for the highest conductivity channel. The 
known resistance vs. output functions determined above were used in conjunction with 
the measurements in standard solutions to determine the probe constant vs. resistance 
relationship The probe constant is calculated by: 
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[ 12 ] 


C P = 


10 6 

Rc 


where Cp is the probe constant [cm], R is the measured resistance [Q], and a is the 

electrical conductivity of the sample [pS/cm] . The probe constant vs. resistance 
functions are shown in Figure 1 1 . Only the lowest conductivity solution (DDI water) fell 
within the useful resistance range of the low a channel, so a cell constant vs. measured 

600 resistance relationship could not be determined for the low a function. The single point 
cell constant has been applied to all low c function data described here. 


E. H 2 0 Vapor 

Calibration of the TECP flight units for humidity was done at the University of 
Washington Mars Atmospheric Simulation Chamber (MASC) facility in March 2006. 
605 Over 50,000 humidity measurements at 5 to 30 second intervals were collected over a 
period of a week. The measurements covered a wide range of frost point temperatures 
(193 K - 263 K) and instrument temperatures (7* 208 K to 303 K). Most of the data was 
collected under 500 [Pa] of CO 2 (7 [mb]), flowing through the chamber at ~200 cm 3 min' 
1 (STP). The absolute humidity was measured by two different NIST-traceable chilled 
610 mirror hygrometers (an EdgeTech DewPrime I and a Buck Research CR-1). 

The air temperature in the chamber was typically a few degrees colder than T h , due to 
self-heating of the circuit board. This temperature difference, combined with the small 
anisothermality of the environmental chamber, prevented us from reaching higher relative 
humidities due to ice formation in the chamber. For example, the vapor pressure 
615 corresponding to R H = 0.60 at -55°C is equal to the saturation vapor pressure at -58.9°C. 
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The ADC counts are related to Tb and the relative humidity, Rh, via a quadratic equation 


ADC = -1.5332 Rl + (963.34 + 14.7057;) R H + (2804.7 - 1.03737; - 6.9089x1 0“ 3 7; 2 ) [13] 

This relationship was obtained by a least-squares fit of nearly 12,000 data points using a 
6 parameter quadratic function of relative humidity and temperature, and has a standard 
620 deviation of 4.5 ADC counts (Figure 12). Applying the calibration involves plugging the 
board temperature into equation 13, and solving for the root of the quadratic. 

The sensor output voltage is measured with a 12-bit ADC, and has a minimum value (at 
Rh= 0) of around 2800 digital number (DN), and an apparently temperature-dependent 
maximum (at Rh= 100%) of at least 3200 DN (the highest value we observed). At 213K 
625 the maximum value of the response function is only 2900, so this standard deviation 
represents a 5% error. However at 233 K the increased range of response reduces the 
standard deviation to about 2%. 

IV. Laboratory and Field Characterization 

Although the TECP meets the measurement requirements as verified in the lab, 
630 practicalities of field and spacecraft operation will affect the accuracy of the various 
measurements on the Martian surface. This section outlines some of those factors, and 
presents preliminary characterization data from the laboratory, spacecraft testbed 
activities, and field testing in the Antarctic Dry Valleys. Additional characterization 
work is underway to further improve science return from TECP. 
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635 A. Characterization 

1. Needle temperature vs. regolith temperature 

TECP can accurately measure the temperature of its individual needles. However, needle 
temperature may not be identical to the regolith or air temperature. If, for example, the 
thermal conductivity of the regolith is extremely low, the TECP needles may be equally 
640 coupled to the regolith and the TECP body. 

To characterize this effect, an experiment was conducted using three different test 
materials: water stabilized with agar, dry silt loam soil, and EPS. In each instance, the 
material was placed in a thermal chamber and brought to a temperature below ambient. 
The sample was then removed from the temperature chamber, placed in an insulating 
645 container at ambient pressure, and the TECP needles were inserted into the sample. The 
needle temperature, 7j, and sample temperature were monitored until they equilibrated, 
and then recorded. This process was repeated at different sample temperatures ranging 
from well below to well above the TECP temperature . 

It is apparent from the data that the error is linearly related to the temperature difference 
650 between the TECP body (7j) and the needles - both measurable quantities. The slope of 
the linear relationship is, to first order, linearly related to the thermal conductivity of the 
medium in which the TECP needles are inserted. Fortunately, thermal properties are 
based on AT measurements, and are independent of sample temperature. So, the 

preliminary characterization data presented here indicate that regolith temperature can be 
655 derived from measured needle temperature, 7j, and the independently measured regolith 
thermal conductivity. 
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2. Imperfect needle insertion 

In an idealized experiment, the TECP needles are fully seated in an infinite homogeneous 
half space, with no voids either at the sides or at the end. However, in practice, the TECP 
660 will be put into use at the end of the 2.3 m Robotic Arm, and the possibility of gaps 
exists. Can measurements be made with imperfect placement? 

A series of experiments were conducted to understand the errors introduced into the 
various TECP measurements from: a) incomplete needle insertion; b) complete insertion 
with partial retraction, and; c) complete insertion with lateral movement. The objectives 
665 were accomplished by using a manual mill to precisely control the position of TECP in 
relation to the sample, and create voids of known dimension. The materials tested during 
the incomplete insertion portion of this exercise were glycerol, stabilized water, AETB-8, 
a thermal protection foam, dry clay soil, dry silt loam soil, and dry quartz sand. All of 
these samples with the exception of the glycerol and stabilized water were tested in the 
670 partial retraction and lateral void exercises. The results of these experiments are 
summarized in Figure 13 for the three soils, and the salient points are descrbed below. 
The description of mitigation approaches against non-ideal insertions is deferred until the 
discussion of TECP-Robotic Arm Interoperability tests (below). 

Electrical Properties 

675 1 . One can think of the overdrive data as an investigation of the effects of compression on 

the sample, the partial retraction as the effects of air gaps, and the lateral movement 
experiment as their combined effects. 
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2. Progressive linear immersion of the needles in the soil particles has a more or less 
linear effect on measured permittivity, but an exponential effect on electrical 

680 conductivity. 

3. Compaction increases the measured permittivity linearly, and the measured 
conductivity exponentially. 

4. The creation of voids around the needles decreases the permittivity linearly, but 
decreases the measured conductivity exponentially. 

685 5. In lateral motions, the effects of compaction initially dominate the electrical 

properties, but beyond 0.25mm of lateral movement, the effects of the voids appears to 
dominate. 

6. The sensitivity of the measured electrical properties to the placement of the needle is 
proportional to their permittivity and conductivity. It’s useful to look at the actual 
690 measured quantities, to see how the effects of insertion gaps are proportional to the 
permittivity and conductivity. 

Thermal Properties 

1. Although the heat capacity measurement is nicely behaved, varying linearly with 
insertion, and being relatively insensitive to compaction, the thermal conductivity results 

695 appear to be complicated by the data reduction technique. 

2. Compression of the soil by the TECP results in linear increases in thermal properties, 
at least up to 3 mm, due to the loss of pore voids in the vicinity of the needles. 

3. Voids result in linear decrease in the measured heat capacity, but the effects on the 
derived thermal conductivity are difficult to deconvolute from the idiosyncricies of the 

700 mathematical modeling that must be done to extract K. 
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4. In lateral motions, the effects of compaction initially dominate both the electrical and 
thermal properties, but beyond 0.25 mm of lateral movement, the effects of the growing 
voids become increasingly important. 


3. Wind 

705 There are no requirements for TECP to measure wind velocity. However, application of 
a protocol similar to the thermal properties experiments has shown some promise in 
characterizing wind speeds, particularly in high-wind regimes where the wind-sock 
associated with the Meteorology package is not capable of measuring wind as well. 

The procedure involves heating needle one for a period up to 70 seconds, and then 

710 following the cooling curve for another 90 seconds. Preliminary analysis of the heating 
and cooling data acquired at Mars-like pressures in a wind tunnel at the University of 
Michigan suggest that both peak temperatures and the cooling time constant 


Tc = 


l_d T 
Tdt 


[14] 


are sensitive to wind velocity. Characterization of the TECP response to wind velocity is 
715 ongoing, and the determination of whether or not to try wind velocity measurements on 
the Martian surface will be made once those analyses are complete. 


B. Antarctic Dry Valley Field Testing 

Field testing of a flight model unit, TECP010, were carried out in Beacon Valley, 
McMurdo Dry Valleys, Antarctica, during the 2007-2008 International Polar Year 
720 deployment. 
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1. Field Site and Data Collection 


The test site in Beacon Valley (77° 50’ 24” S; 160° 39’ 27”E) is at an elevation of 1400 
m. The annual average temperature is -22 °C. Snowmelt is confined to a layer < 5 cm in 
depth, based upon observations of oxidized soil horizons and salt concentrations. Below 
725 this, H 2 O transport is restricted to thin films or vapor. At the time of the field 
observations, there was no snowcover, and the surface was visibly dry. 

The surface materials are silica-rich sands from the surrounding Beacon Sandstone, along 
with clastic fragments of dolerite, which occurs as intrusive sills throughout the Beacon 
Sandstone. The near surface is very rich in salts, predominantly mirabilite 
730 (Na 2 S 04 » 1 OH 2 O), halite, and a number of minor phases. Rocks and boulders are glacial 

erratics with a wide variety of compositions. Additional analyses of the sampled material 
is pending. 

During the austral summer, the sun does not set, so all measurements were conducted 
with the TECP, and the soil into which it was inserted, either in full sunlight, or shadow. 
735 Shadowing resulted from one of three sources. The TECP was placed near the north 
(afternoon sunlit) side of a pyramidal 30 x 19 x 17 cm rock. A Vaisala HM 34 Humidity 
and Temperature Sensor was placed next to it, and used periodically throughout data 
collection for comparison to the TECP humidity data. 

During the “nighttime” hours, (2100h - 0400h), the sun illuminates the soil and face of 
740 the rock opposite the TECP (Figure 14, left). Between 0400h and 0800h, the entire site 
was shadowed by the steep valley walls. The sun emerged from behind the valley wall at 
~0800h, and for several hours, fully illuminated the TECP, the rock face, and the soil 
being sampled. Between 1700 and 173 Oh, the shadow of the Science Tent, from which 
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the TECP was operated, covered first the soil, then the TECP body itself, and finally the 
745 rock. 

Data collection spanned two periods, January 2, 2008 from ~2100h to 0400h, and the 
following day from 1 540h to 1 940h. Raw DN counts for each channel were converted to 
physical units via the functions published for TECP010 in the TECP Calibration Report. 
For the most part, the data set was nominal, and captured expected patterns in thermal 
750 and electrical properties. However, at least one anomalous and interesting event was 
recorded that merits explicit discussion. 

2. Selected Results and Discussion 

Throughout the initial period of data acquisition, the TECP and soil were shadowed, 
although the rock, and the soil opposite the TECP were continuously illuminated. 

755 Needle temperatures decreased steadily until ~0300h. Beyond that, needle and board 
temperatures appear to be approximately stable, suggesting that the power dissipated by 
the TECP electronics is approximately balanced the heat loss of the TECP to the 
environment. 

The most interesting part of the data are the electrical properties. The Bulk Relative 
760 Permittivity (£b) decreases along with temperature until ~0215h, and then begins a rapid 

increase (Figure 14, right).. The total change from one time series to the next is greater 
than the variation within time series, suggesting that the effect is not due to noisy 
readings. The Electrical Conductivity data exhibit a somewhat different, but correlated 
response (not shown). Both the EC Low and EC Medium channels exhibit a sawtooth 
765 pattern; an abrupt drop in a correlates with the minimum in the £b data, after which, o 

and £b increase sharply throughout the remainder of the night. 
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This is a puzzling observation. At the temperature (265 K - 270 K) and the frequencies 
of the TECP measurement (1 kHz for a; 8 MHz for £b), both a and £b are positively 

correlated with temperature in dry, rock-forming minerals (Oheloft, 1981, Buehler et al. 
770 2003; Najipolo et al., 1998). The electrical properties change slope simultaneously, and 

approximately an hour before the cooling rate of the soil changed. 

Since the measurements were acquired in sublimation till above ice-cemented permafrost 
(depth to ground ice at the site was ~25 cm), and the Rh sensor registered a non-zero and 
stable atmospheric relative humidity in spite of ambient cooling, it is reasonable to 
775 suppose that there existed a population of unfrozen water on the surfaces of the soil 
particles. Normally, however, the volumetric unfrozen water content (0) of soil in a 

closed system decreases as the temperature decreases, which would also be inconsistent 
with these observations. 

It is possible however, to increase £b and a in a soil by adding unfrozen H 2 O molecules, 
780 (open system). The correlation between the unfrozen water content and dielectric is a 
well-known phenomenon; there is even a “universal” relationship, which describes the 
relationship between £b and 0 for a remarkably large variety of soil types (Topp et al. 
1980) 

£ b = 3.03 + 9.300 + 146O0 2 -76.70 3 [15] 

785 This equation assumes that a change in 0 is the only variable affecting £b. If applied to 

the data recorded here, it would suggest an increase in unfrozen water of 0.25% between 
0215h and 0345h (from 4.15% to 4.4%). 
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The humidity data are stable, suggesting the atmosphere is not the source of An 
intriguing possibility, which may be consistent with the large d 2 £b/dl 2 at 0215h, is the 

790 arrival at the needles of additional H 2 O molecules, migrating in thin films from the 
continuously warming far side of the rock. While this hypothesis requires further 
analysis, the unexpected relationships between thermal and electrical properties is 
suggestive of the type of events that might prove particularly exciting in data from the 
Martian surface. C. PIT Testing 

795 Since the measured electrical and thermal properties are so sensitive to the details of 
needle insertion, characterization of the interoperability of the Robotic Arm with the 
TECP were conducted in the Payload Interoperability Testbed (PIT) at the University of 
Arizona, to assess the likely scope of gaps associated with non-ideal insertion. 

The original goal was to conduct all TECP insertions along an angle of attack of <0.08 
800 radians, which is the taper angle of the needles. However, this proved impracticable. 
Movement of the RA in the direction of the TECP needles (z-axis) results in bouncing by 
several mm during start and stop transients, with the principle motions directed along the 
x and z axes (see Figure 15 for an example, and the co-ordinate axes referenced in the 
discussion). In the cohesive Baghouse Dust used in the PIT, this motion leaves needle 
805 tracks in the soil, and effects the quality of the TECP insertion. 

To establish an operational baseline for the performance of the RA/TECP operations, we 
first conducted a manual insertion of the TECP into the Baghouse Dust, in a manner as 
near to ideal as possible . Results for various mitigation strategies were compared 
against the results of the ideal insertion to assess possible strategies for mitigating gap 
810 formation. 
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The mitigating approach that proved most successful was to insert the TECP to the best 
of the RA capabilities, and then run the RASP for 1 5 seconds. The RASP is a small drill- 
like sampling device, included in the Robotic Arm Scoop in order to improve the chances 
of collecting very cold, and very hard, subsurface ice. In highly cohesive soils, the high- 
81 5 frequency, small-amplitude vibration caused by the RASP motor, coupled to the needles 

through the RA and TECP body, results in the collapse of any substantial voids 
surrounding the needles. 

Of course, this technique does produce a alteration of the soil particulate structure from 
its original state, but the introduction of any measurement device will have the effect of 
820 disturbing the medium. Continued characterization efforts are underway to 
quantitatively assess the errors introduced by the non-ideal insertions, combined with 
effects of the RASP to collapse the resulting voids. 

D. RAC Images 

One means of characterizing the cleanliness of any given insertion is by imaging the 
825 TECP with the Robotic Arm Camera (RAC). In Figure 16, the RAC has imaged the 
TECP during insertion tests, as a function of the rotation angle of the TECP relative to 
the RA. The insertion angles range from 30° to 77° (nearly orthogonal to the RA). In 
some instances, the entry point of at least one needle can be seen. As the rotational angle 
gets smaller, it is harder to verify insertion and possible gaps with the RAC. Other things 
830 being equal, we prefer to make measurements in an orientation that permits verification 
of TECP placement with the RAC. In any event, the RAC will image any area sampled 
by TECP both prior to insertion and immediately after retraction. 
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V Surface Operations 

Operations Planning for TECP is somewhat unique for the Phoenix mission, because 
835 TECP is an in situ measurement of the regolith properties, which change with the 
environmental variables. In addition, by trenching in the regolith, we will alter the 
physical structure of the regolith, and disturb some of the properties we intend to 
measure. Because we do not yet know the rate at which trenching operations will occur, 
nor the total depth to ground ice, nor the physical configuration of the trench, nor the 
840 availability of undisturbed materials in the workspace accessible to the RA, actual 
scheduling and targeting decisions will have to be made as part of the strategic and 
tactical processes. However, some considerations are described below that should guide 
planning for TECP measurements during the expected mission. 

The physical and chemical activity of water is dependent upon temperature, and 
845 therefore, expected mission planning for characterization of regolith H 2 O activity should 
focus on the local warmest period during operations. 

Landing occurred May 25, 2008 which corresponds to L s =76°, in the late northern spring. 
We can do the arithmetic very crudely and say there are 0.53842° L s /Sol, so the 
approximate mapping of Mission Sol against L s would indicate that summer solstice 
850 would occur around Sol 26 (Figure 17). This will be earlier than the local thermal 
maximum, and subsequent observations should be made approximately every 5° of L s 
until after the local thermal maximum. Therefore, roughly, H 2 O observation campaigns 
would also occur as follows: 

L s 95° = Sol 35 
855 L s 100° = Sol 44 
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L s 105° = Sol 55 


The ideal would be to do an insertion into the soil during the warmest part of the sol, and 
monitor the temperature, electrical conductivity, dielectric permittivity and Rh, but not 
thermal properties (heating during thermal properties measurements will alter, and 
860 perhaps redistribute H 2 O) on or about sols 26, 35, 44 and 55. 

The desired insertion point, or points, cannot be determined in advance, because we do 
not know the configuration of the trench at each date, the specific shadowing associated 
with the lander, the depth to the ice table and so on. However, the H 2 O mobility 
campaigns should be conducted for an hour or two around the thermal maximum of a 
865 target selected during the tactical process in the sols leading up to the thermal campaign 
sols. 

BThe temperature and thermal properties should be measured in at least one material near 
the beginning (L s ~80°), middle (L s ~ 105°) and near the end (L s ~ 125°) of the mission. 
These correspond roughly to sols 10, 50 and 90. The diurnal maximum surface 
870 temperature will change relatively little over the course of the mission, but the nighttime 
minima will vary considerably more (Figure 14). By targeting the beginning, middle, and 
end of the mission, we will be able to characterize the thermal properties of the accessible 
materials, as a function of temperature, over the greatest possible range. 

The temperature and thermal properties measurements should be conducted in material 
875 that is minimally disturbed. If horizontal compositional layering is evident in the dry 
overburden, each layer should be characterized. 
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At least one of these thermal properties campaign days should be correlated with 
THEMIS observations of the Phoenix landing site, to characterize the relationship 
between in situ thermal properties measurements and those derived from remote sensing. 

880 Thus, a reasonable but modest regolith characterization campaign could be conducted by 
prioritizing TECP measurements on approximately 7 sols throughout the mission. The 
specific sols on which the H 2 O and thermal properties are characterized can be decided 
tactically, but the rough guidelines given above will provide good coverage of 
temperature, and maximize the probability of catching any translational or rotational 
885 mobility in regolith H 2 O molecules. 

VI Summary 

The TECP will provide an exciting opportunity to measure the thermal properties of the 
Martian regolith in situ, as well as to characterize the response of the regolith in the high 
latitudes to the extreme variations in boundary conditions, both for temperature and H 2 O 
890 abundance, that are characteristic of the polar summer. In combination with other 
Phoenix instruments, it is hoped that the physical and chemical processes associated with 
energy and H 2 O fluxes through the regolith can be quantitatively understood in the 
present climate, and that this knowledge will give us the capability to decode the geologic 
history of the Martian polar latitudes as uncovered by the Phoenix mission. 

895 TECP is a thoroughly calibrated instrument, but there are practicalities of field and 
spacecraft measurements that must be taken into account. Mitigating procedures, such 
as using the Phoenix Robotic Arm RASP to collapse voids around the TECP needles are 
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described. TECP was field tested in the Antarctic Dry Valleys, and showed promise in 
detecting small changes in unfrozen water under Antarctic conditions. 

900 
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Table 1. Thermal properties standard materials used in ti 

he TECP calibration. 

Material 

Thermal 
conductivity 
[W m 4 K 4 ] 

Source 

Heat 
capacity 
[kJ m’ 3 K 4 ] 

Source 

Macor (glass 

ceramic, 25 C) 

1.402 

Decagon 

2010 

Decagon 

Water (stabilized, 
25 C) 

.608 

CRC 

4174 

CRC 

HDPE (23 C) 

0.513 

Decagon 

1690 

Decagon 

Glycerol (20 C) 

0.285 

Ganic & Hicks 

2963 

Ganic & Hicks 

Expanded 

polystyrene (EPS, 
20 C) 

0.033 

Laser Comp 

63 

calculation 


Table 2. Thermal conductivity (k) and heat capacity (CV) values measured at needle 1 
and 2 with the TECP flight unit. The values of K m and CV, m are generated from equations 
13 and 14, and the known values of K a and Cv, a are taken from literature or other sources. 


Material 

K/i 

192 

C/1 

C/2 

K m 

Cv,m 

Ka 

Cv'd 

EPS 

0.110 

0.199 

1.741 

1.539 

0.033 

0.065 

0.033 

0.063 

EPS 2 

0.112 

0.202 

1.727 

1.545 

0.034 

0.061 

0.033 

0.063 

Glycerol 

0.498 

0.526 

1.815 

5.149 

0.283 

2.913 

0.285 

2.96 

Glycerol 2 

0.517 

0.534 

1.669 

5.097 

0.291 

2.823 

0.285 

2.96 

Water 

0.936 

1.071 

1.077 

7.232 

0.595 

4.296 

0.606 

4.17 

Water 2 

0.951 

1.125 

1.086 

7.348 

0.614 

4.367 

0.606 

4.17 

Water 3 

0.942 

1.035 

1.131 

6.728 

0.602 

4.097 

0.606 

4.17 

Macor 

1.629 

2.846 

0.409 

3.095 

1.391 

2.008 

1.401 

2.02 

Macor 2 

1.636 

2.832 

0.401 

3.149 

1.392 

2.059 

1.401 

2.02 

Macor 3 

1.682 

2.947 

0.351 

3.085 

1.441 

2.024 

1.401 

2.02 

Macor 4 

1.634 

2.904 

0.377 

3.151 

1.404 

1.986 

1.401 

2.02 

HDPE 

0.683 

1.062 

1.268 

3.484 

0.501 

1.679 

0.513 

1.691 

HDPE2 

0.694 

1.122 

1.199 

3.627 

0.514 

1.693 

0.513 

1.691 
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Table 3. The materials and results used in the TECP dielectric calibration 


Material 

Nominal 

dielectric 

Temperature 

(C) 

NIST Dielectric 

Output TECP 
FM (ADC) 

Air 

1 

23 

1 

642 

1 -decanol 

7.5 

24 

8.07 

1410 

Tert butanol 

10.9 

32 

11.22 

2300 

Isopropyl 

18.3 

23 

19.76 

4016 


Table 4. Useful resistance ranges for each TECP measurement range 


Channel 

minimum resistance 

maximum resistance 

Hi Conductivity (ECH) 

— 

6kQ 

Medium Conductivity (ECM) 

6kO 

1MQ 

Low Conductivity (ECL) 

156kO 

-- 
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Table 5. The coefficients for the 4 th order polynomial fit of ln(7?) to the ADC counts as a 
function of temperature. 


T[ K] 

Hi range 

160 

In ( R ) = -8.341x 10' I4 ADC 4 + 8.237x10'° ADC 3 - 2.75 lxlO 06 ADC 2 + 4.743xlO-° 3 ADC + 3.136 

200 

In ( R ) = -7.892x1 O’ 14 ADC 4 + 7.927x1 0'° ADC 3 - 2.691x1 O' 06 ADC 2 + 4.71 8x1 O' 03 ADC + 3.128 

240 

In ( R ) = -7.809x1 O' 14 ADC 4 + 7.858x1 0 10 ADC 3 - 2.673x1 0 06 ADC 2 + 4.699x1 O' 03 ADC + 3.136 

280 

In ( R ) = -7.742x1 0 14 ADC 4 + 7.807x10'° ADC 3 - 2.66 lxl O' 06 ADC 2 + 4.690x1 O’ 03 ADC + 3.139 

323 

In (i R ) = -7.769x 10" I4 ADC 4 + 7.820x1 0 10 ADC 3 - 2.662x1 O' 06 ADC 2 + 4.689x1 0 03 ADC + 3.141 


Mid range 

160 

In ( R ) = -3.248x1 0 14 ADC 4 + 4.654x1 O' 10 ADC 3 - 1.922x1 O' 06 ADC 2 + 4.037x1 O' 03 ADC + 7.918 

200 

In (R) = -3.248x1 O' 14 ADC 4 + 4.654x1 0 10 ADC 3 - 1.922x1 O' 06 ADC 2 + 4.037xl0-° 3 ADC + 7.918 

240 

In (R) = -3.263x1 0 14 ADC 4 + 4.605x1 O' 10 ADC 3 - 1.903x1 O' 06 ADC 2 + 4.019xl0-° 3 ADC + 7.905 

280 

In (R) = -3.128x 10 14 ADC 4 + 4.498x1 0 10 ADC 3 - 1.877xl0-° 6 ADC 2 +4.000xl0° 3 ADC + 7.901 

323 

ln(fl) = -3.100xl0 14 ADC 4 + 4.478x1 O' 10 ADC 3 - 1.873xl0-° 6 ADC 2 +3.998xl0 03 ADC + 7.901 


Low range 

200 

In (R)= 1.3236x10 s ADC 3 - 1.1648xl0 4 ADC 2 +3.4326x10 ' ADC + 323.41 

240 

In (R)= 1. 21 65x1 0 8 ADC 3 - 1.0593xl0 4 ADC 2 +3.0932x10 ' ADC + 286.57 

280 

In (R)= 1.9291x 10 8 ADC 3 - 1. 7078xl0 4 ADC 2 + 5.0553x10 ' ADC + 483.93 
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Figure 1 



Phoenix Region A-D, Box 1 
Geomorphic Map Units 

H Highland 

H Lowland - Dark 


Lowland - Bright 
Knobby Terrain 
Intermediate Terrace 


Blocks/Mesas 


Crater Interior 


Crater Ejecta 

Map colors are overlain on a 
THEMIS visible mosaic, with MOLA 
shaded relief providing extended 

K.D. Seelos 
JHU/APL 
March 4, 2007 


Figure 1 . Geologic Map of the Phoenix Landing site, by K. D. Seelos (see Arvidson et al., 
This Issue), with examples of patterned ground seen in Hi-Rise imagery (all examples ~ 50 
cm/pixel). 
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Figure 2 



Humidity Sensor 





5 \\ 


Reference T for thermal measurements 
Ground for dielectric measurement 

Dielectric measurement excitation 


AT for C and k measurements 
EC measurement excitation 


Heated for thermal measurements 
AT for thermal conductivity measurement 
Ground for EC measurement 


Figure 2. Photograph of the TECP instrument (top) and with the external cover 
removed to allow access to the electronics board (bottom). For each needle, the 
numerical designation and functionality is identified at right. The TECP is 1 18.76 
mm in length. Its needles are 1 5 mm long. 
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Figure 3 


Electrical Conductivity Dielectric Permittivity 



Temperature/Thermal Properties 


Humidity 


thermocouples 

(AT) 



Figure 3. The electrical circuits employed by TECP to measure electrical and thermal 
properties of the regolith and atmosphere. 
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Figure 4 


1075 


1080 



Figure 4. The TECP mounted on the Robotic Arm. The Robotic Arm Camera (RAC) is 
at right, and can image either the scoop contents, or the TECP when it is rotated counter- 
clockwise for use. 
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Figure 5 



Board ADC 

Figure 5. TECP FM board temperature sensor calibration. 
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Figure 6 



Figure 6. The calibration setup for the thermocouple calibration. 
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Figure 7 



T (K) 


Figure 7. The Seebeck Coefficient of the thermocouple wire used for the TECP FM. 
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Figure 8 



Cv.mfMJ/rn^K' 1 ) 


Figure 8. The volumetric heat capacity calibration curve. C v , a is the volumetric heat 
capacity of the laboratory standard materials, and C v>m is the measured heat capacity - 
which is actually derived not through direct measurement, but through fitting Eqs. 10 - 
1 1 to the measured heating and cooling curves associated with heating one of the needles. 
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1095 Figure 9 



Figure 9. The thermal conductivity calibration curve. Ka is the thermal conductivity of 
the laboratory standard materials, and K m is the measured thermal conductivity - which is 
actually derived not through direct measurement, but through fitting Eqs. 10-11 to the 
1100 measured heating and cooling curves associated with heating one of the needles. 
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Figure 10 



ADC 

Figure 10 . The dielectric permittivity calibration curve for the TECP flight unit. 
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Probe Constant [cm] 


1105 Figure 11 



In R 

Figure 11 . The Probe Constant for the TECP as a function of the resistance, for each 
of the 3 channels of the EC function. 
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Figure 12 


Humidity Calibration Function 



Humidity Calibration Data Points 



Figure 12. (A) Contours of the calibration function for the TECP humidity sensor at four 
different values of T b . The dashed portions of each line indicate an extrapolation beyond 
the range of data used to derive the calibration function. (B) Scatterplot of the 
environmental conditions corresponding to each calibration data point. 
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Figure 13 
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Figure 14 




Time (D:HH:MM:SS) 


Figure 14. The TECP deployed in Beacon Valley (left). Needle 3 temperature and bulk 
dielectric permittivity recorded during the overnight hours of 2 January, 2008. 
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Dielectric Permittivity 


Figure 15 



Figure 15. The TECP during interoperability testing with the Robotic Arm in the 
Payload Interoperability testbed. The axes refer to the TECP-centered co-ordinate 
system used by the RA. The needle holes are evidence of overshooting in the z 
direction during a commaned move of + 1 cm in the z direction. Their shape reflects 
uncommanded movement of ~3.5 mm along the x axis. 
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Figure 16 



Figure -16 . The RAC view of the TECP at insertion angles of 30°, 47° and 77° (top 
to bottom'). 
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1115 Figure 17 



130 160 190 220 250 


Figure 17. The predicted surface temperature at 70°N (top), The arrow and vertical 
red lines indicate the nominal Phoenix mission duration. For comparison the lower 
plot is Tig, indicating approximate surface kinetic temperature. (From Kieffer and 
Titus, 2001) from TES. The horizontal white line indicates 70°N. 
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